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Abstract 

We have calculated the radiative corrections to the decay r — ► n(K)v T , taking to ac- 
count internal bremsstrahlung and structure dependent radiation in the radiative decay 
and point meson, hadronic structure dependent and short distance contributions in the 
virtual corrections. We display the spectra of the photon energy and of the pion-photon 
invariant mass in the decay r — ► iri/ T j and compare with the PHOTOS Monte-Carlo. 
Our result for the radiative correction to the ratio T(t — ► 7r^ r (7))/r(7r — ► ^^(7)) is 
SR t/7T = (0.16 ± 0.14) % and for the ratio T(r -> Ku T (j))/T(K -> ^(7)), we obtain 
SR t/k = (0.90 ± 0.22) %. 
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1 Introduction 

1.1 Lowest Order Prediction 

The simplest semihadronic decay mode of the tau lepton is unique in that there is essentially 
no uncertainty in the theoretical prediction at the order 0(a°). The decay width T(t — » tiu t ) 
can be predicted using the precisely measured pion decay width r(7r — > \iv^). The only 
parameter in both decays is the pion decay constant f n defined by the matrix element of the 
weak hadronic current, 

<0\J^ ak \7T+(p)>=iV2f^ (1) 

In the ratio Rt/ w of the decay rates 

r(7T -> fiU„) 

the coupling /„. cancels, and the resulting formula involves the masses of the particles only, 

1 m 3 T (1 - ml/ml) 2 

2 m^mt (1 — mrjm%) 

Multiplying this ratio by the experimental decay width r(7r — ^ A"^)> one § e ts the 0(a°) 
prediction for Y(r — ► 7rz/ r ). 

1.2 Radiative Corrections of the order 0(a) 

Because of the infra-red divergences, at the order 0(a) one has to consider the ratio of 
inclusive decay rates into final states with or without an additional photon: 

R _ T(t -> 7rz/ T ( 7 )) _ 1 ml (1 - ml/m 2 T f , | §R x 

r(7T — >• HVfid)) 2m n m 2 l (l — m^/rn^) 2 ^ T ' 

The radiative correction 5R T / n arises from the general Feynman diagrams shown in Fig. [I], 
and the corresponding diagrams for the decay n — ► fii/^j) . There are the virtual corrections, 
Fig. [I] (a)-(c), and the corrections due to the radiative decay, (d)-(e). The diagrams fall into 
three classes. The diagrams of the first class (Fig. [I] (a) and (d)), where the photon only 
couples to the leptonic side, are determined by f n and QED, and so they can be calculated 
without theoretical uncertainty. The diagrams of the third class (Fig. [I] (c)), where the photon 
couples twice to the hadronic side, give identical corrections for the tau and pion decay and 
so cancel in R T / n . Thus the essential problem is to calculate the diagrams of the second class, 
where the photon couples once to the hadronic side, Fig. [I] (b) and (e). 

What makes the calculation of these radiative corrections difficult and interesting is the 
fact that due to the integration over the momentum of the virtual photon, the physics of very 
different energy scales Q is involved. For very small energy scales, Q 2 -C m 2 the relevant 
amplitudes are fixed by low energy theorems of QCD. For intermediate scales, Q 2 ~ m 2 the 
amplitudes are dominated by hadronic resonances, and for large scales Q 2 3> m 2 , they are 
dominated by the short distance behaviour of the weak interaction. 

In Sec. 2 we will discuss the radiative decay, which is not only an essential part of the 
total radiative correction, but also of interest in its own right. In Sec. 3 we will then consider 
the virtual corrections and the total correction to the ratio R T / n . 



Figure 1: Diagrams determining the radiative corrections to r — > ixv T 
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2 The Radiative Decay r — ► 7rz/ r 7 

The amplitude for the radiative decay r — ► 7rz/ T 7 can be divided into the so-called internal 
bremsstrahlung (IB) and the structure dependent radiation (SD) |I|, |2|, |3|, [|, [^]: 



M\t-(s) 



it {p)v T {q)l{k)} = M 



IB 



M 



SD 



(5) 



The internal bremsstrahlung corresponds to hooking photons to the external lines, considering 
the pion as an structureless elementary particle: 



Mib = -G F cos 6 c ef n m T u u (q)^ 



p ■ e 
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s ■ e 



p ■ k 2s ■ k s ■ k 



u T (s) 



(6) 



The structure dependent radiation can be defined as the contribution from the diagram in 
Fig. Me) minus the effective point pion contribution, resulting in 



M SD = -^^(W^7-^W^ 



V2 



m„ 






m T . 



(7) 



This amplitude M.sd involves two form factors Fy(t) and F^(t), which depend on the invariant 
mass squared of the pion-photon system, t = (p+k) 2 . The vector and axial vector form factors 
F v and Fa correspond to the J p = 1~ and J p = 1 + projections of the W boson, respectively. 
And so in order to be able to calculate the decay rate, we have to parameterize the form 
factors F v (t) and F^it). 



2.1 Normalization at t — 

We use the following values for the form factors at zero momentum transfer: 

F A (t = 0) =0.0116 ±0.0016 (8) 

Here Fy(Q) is obtained from the axial anomaly, whereas F A (0) has been measured in the 
radiative pion decay n — > ev e ^ ||. 

However, for the radiative tau decay, we need Fy and Fa in the whole range m 2 < t < m? T . 
In order to parameterize the energy dependence of the form factors we assume vector and 
axial vector meson dominance. 

2.2 Parametrization of F^it) 

In the case of the axial form factor, the only established resonance which has the correct 
quantum numbers is the ai(1260). Thus we assume 

F A (t)=BW ai (t)F A (0) (9) 

where in the normalized Breit Wigner resonance factor 

2 
777 

m\ - t - im ai V ai (t) 

we use an energy dependent width T ai (£) as calculated from the a± —* 3n and a± — > pn phase 
space 0. 

2.3 Parametrization of Fy 

In the case of the vector form factor, there are three resonances with the correct quantum 
numbers in the mass range covered by the tau, viz. the p = p(770), p' = p(1450) and the 
p" = p(1700). And so the vector meson dominance ansatz for Fy(t) reads 

F v(t) = 1 fivr f BW pW + ABWp/(t) + pBW p „(t)} (11) 

1 + A + p L J 

We have to fix the two parameters p and A which determined the size of the contribution of 
the higher radial resonances. Noting that Fy(t) is related by CVC (conserved vector current) 
to the form factor F nll ±(t) which determines the coupling 7r°77*(g 2 — t) ||, we are able to 
find four constraints for the two parameters: 

1. The limit of F 7r77 * for large negative t is predicted by perturbative QCD || 

hmjF 7Tir {t) = 2f 7T (12) 

2. The slope of F 7r77 *(t) at t = has been measured in 7r° — > e + e~7 0. 



1 7" — >Z/7T ^"^ 



Figure 2: Photon energy spectrum of the decay r — > nUrj: Total prediction (solid), pure 
internal bremsstrahlung (dashed) and the prediction from TAUOLA + PHOTOS (dots with 
statistical error bars) 
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3. The value Fy(m 2 ) of the vector form factor at the rho mass is related to the decay width 
T(p —>■ 7T7) which has been measured. 

4. Similarly, Fy(m 2 ,) is related to T(p' — > nj), which has not been measured but is in turn 
related by vector meson dominance to the measured width T(p' — > ixu). 

The parameter choice 

A = 0.136; fi = -0.051 (13) 

simultaneously fulfills all four constraints. 

2.4 Numerical Results 

Having obtained parameterizations for the hadronic structure dependent form factors, we 
can calculate decay distributions for r — > 7rz/ T 7 Q. In Fig. |2| we display the photon energy 
spectrum dT/dx of the decay. The dimensionless parameter x is defined by 

2.£L 

x = — 1 14 

m T 

where E 1 is the photon energy in the tau rest frame. In Fig. |2|, we do not only display the 
total spectrum, which we predict, but also the contribution from pure internal bremsstrahlung 
(IB) (i.e. neglecting hadronic structure dependent effects) and the result from the Monte- 
Carlo program PHOTOS |10| . PHOTOS provides a semiclassical approximation of the IB 



contribution. We find from Fig. |2] that PHOTOS gives a very good approximation to the IB 
part. Furthermore we find that the IB is strongly dominating the total spectrum, except for 
the region of very hard photons where the rate is extremely small anyway. 

So in fact we find that PHOTOS gives a very good approximation to the photon spectrum 
of t — > 7rz/ T 7, and that hadronic structure dependent effects are almost invisible in this 
spectrum. 

A observable much better suited to separate the structure dependent effects from the in- 
ternal bremsstrahlung is the pion-photon invariant mass spectrum dT/dm nl , which we display 
in Fig. ??. The dimensionless parameter z is defined by 

z = —l (15) 

In this spectrum, a clear p resonance peak and a smaller Oi peak are visible. 



Figure 3: Pion-photon invariant mass spectrum of the decay r — > -kv t ^\ Internal bremsstrah- 
lung (IB, dashed), structure dependent radiation (SD, dash-dotted), interference of IB and 
SD (INT, dotted), and the total spectrum (solid) 
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A measurement of this spectrum would be very interesting both in order to check our model 
and in order to have another channel to measure the poorly known width T ai . However, there 
is a very large background from the decay r — > 7r _ 7r°z/ r , where the 7r° decays into two photons 
and one of these may escape detection. 

It is interesting to mention the total size of the hadronic structure dependent contribution 
(integrated over the full phase space), 

T S b + int(t ~ Trzyy) = ((m _ QM)% = Qm% (16) 



r (r -> -ku t ) 

where the first number (0.11) corresponds to the structure dependent amplitude squared 
(SD) and the second one (-0.06) to the interference between structure dependent radiation 
and internal bremsstrahlung (INT). 

3 Virtual Corrections 

Calculating the virtual corrections, we have to perform loop integrations over the Euclidean 
photon momentum k\ from k\ — up to k\ = m|. In order to be able to perform this 



integration, we separate the integration range into two parts [[LI], [12], [n|. In the long distance 
region of k% — 0- • ■ ^ 2 ut ~ 0(1 GeV 2 ), mesons are the relevant degrees of freedom. For this 
region we build a realistic phenomenological model, taking into account low energy theorems 
of QCD, assuming vector meson dominance and considering experimental constraints such 
as the electromagnetic form factor of the pion, which is known quite well in the relevant 
momentum region. For the short distance integration, k\ = n 2 cut - ■ -m 2 ?, we use the parton 
model. 

3.1 Long Distance Part 

In the integration over small photon momenta, k% = • • • n 2 cut) we proceed as follows. We start 
from the diagrams which correspond to the lowest order, 0(P 2 ), of low energy QCD (effective 



point pion [|1J], see Fig. [|). We then modify these diagrams to account for the momentum 
dependence of the electromagnetic form factor of the pion, which is known to be dominated 
by rho-like resonances, see Fig. [5] We then add to this the loop diagrams which correspond 
to the structure dependent (SD) part in the radiative decays, see Fig. ||. As indicated in the 
figure, we assume a double vector meson dominance of the relevant form factors Hy and H^, 

H v (k,p) =BW UJ (k 2 )F v [(k + p) 2 } 



Figure 4: Effective Point pion diagrams 
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(a) SMi 




(b) 8M 2 
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(c) SM 3 




(d) SM 4 
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(e) SM 5 









Figure 5: Vector meson dominance of coupling of the photon to the pion. In fact we include 
the p, the p' and the p" with relative strengths which have been obtained by fitting the 
measured pion electromagnetic form factor 0. 
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H A (k,p) = BW p (k 2 )F A [(k + pf 



(17) 



where k is the momentum of the virtual photon, p the pion momentum, and Fy and F A are 
the form factors involved in the radiative decays (real photons). In order to estimate the 
model dependence, we also compare to a form with a single vector meson dominance, 

Hy(k, P ) =F V [(k+ P ) 2 ] 

H A (k,p) =F A [{k+pf] (18) 



Figure 6: Hadronic structure dependent loop diagrams 
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Figure 7: Short distance diagrams 
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(a) SA a 




(b) SA b 
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(c) SA C 






(d) SA d 
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(e) SA e 




(f) SA f 







3.2 Short Distance Corrections 

For large virtual photon momenta, k 2 E = fi 2 ut ■ ■ ■ m 2 z , we calculate the corrections to the 
elementary vertex r — > v T ud, see Fig. |7|, and plug the result into Fig. ||. Doing the same thing 
for the pion decay, we obtain the short distance contribution to the radiative correction to 
the ratio R r / n in the form 



6R 



t/tt 



*(l 



3 

2/; 



+i 



du § n (u)r(u) 



(19) 



The integration extends over the scaled relative momentum of the two quarks in the infinite 
momentum frame. $tt(u) is an unknown parton distribution function, and r(u) has been 
calculated from the short distance diagrams. r(u) varies only very little over the integration 
range, and so we can approximate it by its value at u — 0, 



K/*L^(0); 



+i 



du ^-rrlu) 



r(0) 



(20) 



where last equation follows from a sum rule ||15||. In fact we find that the short distance 



correction is dominated by a leading logarithmic contribution |T3|, [16 

2a ml 

' sd 



[SR T /n 



n m? T - n 2 cut 
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f^cut 



m T 



(21) 



3.3 Numerical Results 

Adding up long and short distance corrections, we obtain the total radiative correction, which 
depends on the choice of the matching scale /i cui and on the choice of the parameters of the 



8 



Figure 8: Decay r — > ttv t via an intermediate quark-antiquark state. The bubble on the left 
hand side is to be replaced by the short distance diagrams 5Ai of Fig. 
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hadronics. In Fig. [3.3| we display the radiative correction to the width r —* ixv T in variation 
with ficut, using three different sets for the hadronic parameters. There is a standard set 
(I) with some central values, and sets (II) and (III) with reasonable variations around these 
central values. While for the sets (I) and (III) the dependence of the correction of the matching 
scale fi cut is reasonably small, it is unacceptably large for (II). In fact ,(I) and (III) use the 
double vector meson dominance parameterization of fll~7|), (II), however, uses the single pole 
form of flT8|), which therefore can be excluded. 

We obtain similar results for the radiative correction to r(7r — > {iv^). Taking the ratio, we 
find 

SR t/it = (0.16 ± 0.08 ± 0.04 ± 0.07)% = (0.16 ± 0.14)% (22) 

The first error quoted (0.08%) arises from the matching uncertainties in T(r — > 7rz/ T (7), the 
second one from the matching in r(7r — > fii> ll ('-f), and the last one from the uncertainties in 
the hadronic parameters. The number given here for 5R T / n , (0.16 ± 0.14)%, differs slightly 
from the number we quoted in IJ13J , SRt/k = [0-lG^o!u) %■ This is due to the fact that in 



||13| , we estimated the matching uncertainty by considering the variation of the ratio R T / W , 
whereas now we have discussed the matching uncertainties of the individual decay rates. The 
latter approach should give a better estimate of the true model dependence, because in the 
ratio scale dependences associated with a mismatch of long and short distances tend to cancel 
at large scales. 

Note that the rather small total radiative correction arises from cancellations of larger 
numbers with opposite signs. 

Similarly, we obtain the radiative correction to the ratio of the tau decay into a kaon 
r(r — > Kv T {^)) and the decay width of the kaon T(K — > fiv^) 

SR T/K = (0.90 ± 0.08 ± 0.09 ± 0.14)% = (0.90 ± 0.22)% (23) 

Our result for R T / n should be compared with the recent estimate by Marciano and Sirlin 
17 1, which in terms of R t /-k reads 



SR T/7T = (0.67 ±l.)% (24) 

where the ±1.% is the authors' estimate of the missing long distance corrections to the tau 
decay rate, which they did not calculate. And so we confirm their result within their esti- 
mated error bars, but we are able to reduce the error substantially by performing a complete 
calculation. Translating the radiative corrections into predictions for the branching ratios, we 



Figure 9: Radiative correction to T(r — > ttu t ), using different choices for the parameters of 
the hadronic structure dependent correction: Standard choice (I) (solid), and variations (II) 
and (III) (dashed and dotted, respectively) 
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Hcut [GeV] 



get 



SD + INT 



BR(r -»■ m/ T (i)) =(10.946 + 0.005 ±0.020)% x 

SD + INT 

BR(r -> Ku T (i)) =(0.723 + b~002 ±0.004)% x 



( T r 



V291.6fs 
V291.6fs 



BR(r -► hv T (i)) =(11.669 + 0.007 ±0.021)% x 



r T 



(25) 



.291.6 fs, 

where h denotes the inclusive sum of pions and kaons. Note that here we have separated out 
the contribution SD ± INT associated with structure dependent radiation, corresponding to 
the decay chains r — > pv T , p — > n-f and r — > a\v T , a,\ — > n-f. Whereas we have included all 
photons in 5R T / n , experimental numbers for the branching ratios do not include hard photons. 
Note, however, that the size of this SD + INT part is in all three cases small compared to the 
overall uncertainty of the prediction. Using the new world average for the tau lifetime [|18] 



r T = (291.6 ± 1.6) fs 



(26) 



we obtain 



BR 



theo 



(10.95 ±0.06)% 
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BR^ eo = (0.723 ±0.006)% 



o 



BR™ eo = (11.67 ±0.06)% (27) 

These predictions agree within one standard deviation with the new world averages [|1|] quoted 
at this conference. 
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